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Abstract 

Adding separate chemical potentials A and 7 for Z2-monopoles and vor- 
tices respectively in the Villain form of the mixed fundamental-adjoint action 
for the SU(2) lattice gauge theory, we investigate their role in the interplay 
between the deconfinement and bulk phase transitions using Monte Carlo 
techniques. Setting A to be nonzero, we find that the line of deconfinement 
transitions is shifted in the coupling plane but it behaves curiously also like 
the bulk transition line for large enough adjoint coupling, as for A = 0. In a 
narrow range of couplings, however, we find separate deconfinement and bulk 
phase transitions on the same lattice for nonzero and large A, suggesting the 
two to be indeed coincident in the region where a first order deconfinement 
phase transition is seen. In the limit of large A and 7, we obtain only lines of 
second order deconfinement phase transitions, as expected from universality. 
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I. INTRODUCTION 



The continuum limit of a lattice regularized gauge theory is defined at its critical point 
where the lattice correlation length is infinite. One therefore expects that a large number 
of lattice actions, differing from each other by only irrelevant terms in the sense of the 
renormalization group, describe the same continuum physics. 

For pure gauge theories, the simplest discretized action, the Wilson action |l||, is widely 
used in both analytical and numerical studies. It has been very successful in revealing 
several nonperturbative features of gauge theories. Most such investigations are, however, 
necessarily carried out for a finite value of the lattice spacing a. It seems therefore imperative 
that the universality of these results is verified by employing other forms of lattice actions. 
Such a study of the universality of the deconfinement transition for SU(2) gauge theory was 
initiated in Ref. for the Bhanot-Creutz action || 

Sbc = P f E (l - \^fU P ) + Pa E (l - \Tr a U P ) , (1) 

where the summation runs over all the plaquettes of the lattice and the subscript, f(a) 
indicates that the trace is taken in the fundamental (adjoint) representation. The Wilson 
action corresponds to setting j3 a = in Eq. ([[]). It has a second order deconfinement phase 
transition at which the order parameter (\Lf\) acquires a nonzero value, where 

1 N T 

L f (x) = ~Tr f 1[U t (x,t) , (2) 

T = l 

and U T are the gauge fields in the time direction. The finite size scaling analysis |4[] of 
its susceptibility, X\L f \ = (£/) _ (\Lf\) 2 , yielded an exponent u = 1.93 ± 0.03, in good 
agreement with the corresponding value (1.965±0.005) for the three dimensional Ising model. 
Universality of the continuum limit predicts a similar deconfinement transition belonging to 
the same universality class as the three dimensional Ising model for all values of j3 a . 

Monte Carlo simulations for Eq. (H) on N% x N T lattices with N T = 4 showed || 
that while the second order deconfinement transition point for the f3 a = Wilson action 
entered the f3j — (3 a plane as a line of second order transitions, the transition, surprisingly, 
turned first order for large enough f3 a - This was evident from the facts that i) the order 
parameter, became nonzero discontinuously at the transition and ii) the exponent 

of its susceptibility changed from the Ising model value to 3. If the change of the order of 
the deconfinement transition were to persist in the continuum limit, it would be a serious 
violation of universality. On the other hand, the deconfinement line was also found to 
coincide with the known bulk transition line || for N T = 4. Studies with varying N T further 
revealed that it scarcely moves as one increases N T , especially in the region where a strong 
first order deconfinement transition is observed |J. 

Similar results were also obtained in studies with a Villain form of the mixed 
fundamental-adjoint action [0, 

S = PfE( 1 - l TT f U p) + & E (l - ^P-Tr/t/p) , (3) 
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where a p are auxiliary Z 2 variables defined on the plaquettes and the partition function has 
a sum over all possible values of the a p variables. This summation ensures that the second 
term gets contributions only from the integer representations of SU(2). The deconfinement 
phase transition in the f3f - (3 V plane of this action was shown || to have similar features as 
above : the deconfmement transition turned first order for large enough (3 V and the transition 
line merged with the known bulk transition line for this action j7|. 

Another example of similar interplay of bulk and deconfmement effects was shown []5],|9| 
to be the SO (3) lattice gauge theory, i.e, for f3f =0 in Eqs. (p]) and fl3|). Being in the 
same universality class, it should also have a second order deconfmement phase transition, 
although it does not have an order parameter like the SU(2) theory. The only transition 
found in Refs. ||J| for different N T coincided with the known first order bulk transition. 

One way the above conflicting indications can be interpreted is to postulate that the 
bulk transition interferes with the deconfmement line in such a manner as to make the two 
lines coincident for a substantial range of N T and couplings for all the above actions. A 
possible way to test this hypothesis is to modify these actions further by such additional 
irrelevant term(s) that the bulk transitions are either shifted in the coupling plane or en- 
tirely eliminated. The bulk transition in the Villain form of the SO (3) gauge theory is 
known to be caused by a condensation of Z 2 monopoles in the strong coupling phase ||TU| . 
These monopoles are absent in the weak coupling region, and can be suppressed at stronger 
couplings by adding an irrelevant term to the action. Investigations [|ll|] of the monopole 
suppressed theory at finite temperature found a continuous deconfmement transition, very 
similar to the transition for the SU(2) theory with Wilson action. This suggests trying an 
extension of the technique of suppression of bulk transitions to the mixed actions as well. 
For the action in Eq. (Q), the bulk transitions are caused by the Z 2 monopoles and vortices 
[7J. In this work, we study its finite temperature phase diagram with suppression of these 
objects. 

The plan of our paper is as follows. In the next section, we review briefly the bulk 
phase diagram for the action (|3]) and the results of studies of deconfinement transition for 
it. Following the SO(3) results, we suppressed only the Z 2 monopoles for it and studied 
the resulting theory in the (3f - (3 V plane. As seen in Sec. [TTT], this resulted in suppression 
of the bulk transition line coming from the SO (3) axis, but it was found inadequate for 
suppressing the remaining bulk transition line, where again the same interplay of bulk and 
deconfinement effects was seen. Nevertheless, we were able to demonstrate that a second 
order deconfinement phase transition and a first order bulk transition both exist on the same 
lattice but at different locations in a small range of /3 V . In Sec. |V[ we added an additional 
term to suppress the Z 2 vortices. This resulted in a complete elimination of all the bulk lines 
and a universal deconfinement transition in the entire (3f - f3 v plane was found. A summary 
of our results and our conclusions are presented in Sec. [V]. 



II. PHASE DIAGRAM FOR MIXED VILLAIN ACTION 

The phase diagram of the mixed Villain action, Eq. (||), was studied in Ref. |7j] on 
symmetric lattices and was found to be qualitatively similar to that for the Bhanot-Creutz 
action, Eq. ([]]). As shown in Fig. [I], lines of first order transition, emanating from the 
j3 v — > 00 region and the /3f — axis, join, and extend to the low j3 v phase, before ending 
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at f3 v ~ 2.2. The bulk transitions in this action can be understood to be due to certain Z 2 
topological objects. Defining the monopole density, M, and the electric current density, E, 
as 

M = 1 - ( -jjr $>c ) with a c = J] a p 

\ iV c c / p& dc 

and E = 1 - ( tj E °e / witn °"e = II °p ( 4 ) 



where iV c and iVj are the number of elementary cubes and links of the lattice, respectively, 
the bulk phases are characterized by presence of condensates of these objects J7J. The 
strong coupling phase has both E and M nonzero, corresponding to condensation of both 
Z 2 electric vortices and magnetic monopoles, while the large f3f and (3 V phase, relevant for 
the continuum limit, is free of both condensates, as summarized in Fig. |l] (~ in the figure 
means zero up to 0(exp(— (3f)) ). 

The finite temperature phase transitions for this theory have been explored in Ref. pi. 
As mentioned in Sec. |, the results are similar to those obtained in Refs. 0J^] for the Bhanot- 
Creutz action. The second order finite temperature transition line, coming out of the f3 v = 
axis, turns first order and joins the bulk transition line at its endpoint fl2fl . The finite 



temperature transition on the (3 V axis was studied in Ref. M . It was found to be a coincident 
bulk and deconfinement transition for N T = 4, 6 and 8 lattices. Since the bulk transition on 
the /3 V axis is due to the melting of the monopole condensate present in the strong coupling 
phase (Fig. |l|), and since the monopoles are not present in continuum limit, we studied the 



deconfinement transition on this axis after suppressing the monopoles pT| . The resulting 
theory i) had no bulk transitions, ii) yielded a deconfinement transition very similar to that 
found for the Wilson action, in agreement with expectations based on universality and iii) 
displayed a shift in the transition point with N T . As a natural next step, we turn to a study 
of the mixed theory (action @) with monopole suppression, as described below. 



III. PHASE DIAGRAM WITH MONOPOLE SUPPRESSION 



Adding a chemical potential term for the monopoles as in Ref. |TT |, the mixed 
fundamental-adjoint Villain action becomes 

S(U,a) = P f Y,(l-\ r Tr f U p ) (5) 



v 



where the last summation runs over all the elementary 3-cubes of the lattice, and a c is defined 
in Eq. (^). The classical continuum limit is left unaffected by the additional monopole term 
and one obtains the same continuum relation, 4g~ 2 — /?/ + (3 V , in that limit. For SO(3), A = 
1 was found to be sufficient to suppress the bulk transition. We therefore took A = 1 for our 
simulations in the entire /3f - (3 V plane. Our simulations consisted of a 3-steps iteration. First, 
all the gauge variables were updated using Creutz's heatbath algorithm. This was followed 
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by a heatbath sweep for the Z 2 variables. In the third step of our iteration, a fraction of the 
links (arbitrarily chosen to be 7) were multiplied by a Z 2 element subject to a probability 
determined by the j3f term. This third step is essential for reducing the otherwise enormous 



autocorrelations for large A simulations, and generalizes the similar step used in Ref. JTT 
for the SO (3) gauge theory. Measurements were made after every such compound iteration. 
Using hysteresis runs of 15000 iterations per point, and monitoring the plaquette variables 
P = (\TifUp) and P a = (^a p TrfU p ), we mapped out the phase diagram on an 8 3 x 4 lattice. 

First order transitions, with significant discontinuities in the average plaquette P for 
/3 V > (3f and the "adjoint" plaquette P a for (3 V < (3f, were observed. The transition points 
are shown by filled circles in Fig. 0. While the transition line in the large (3 V region is 
very similar to the A = case (Fig. [I]), unlike that case the transition line does not have 
an endpoint and divides the (3f — (3 V plane in two disjoint parts. Also the transition line 
coming out of the (3f = axis in Fig. [l] is absent here. Since it is caused by condensation 
of monopoles, which have been suppressed here, its absence was to be expected. 

To look for the deconfinement transition, we monitored the behavior of (\Lf\). The de- 
confinement transition point on the Wilson axis was seen to extend into a line of continuous 
deconfinement transitions as we switched on f3 v . On increasing ft v , the deconfinement tran- 
sition line was seen to merge with the first order line, and the order parameter (\Lf\) showed 
a discontinuous jump to a nonzero value, indicating a first order deconfinement transition 
for these points. The dotted (solid) line in Fig. || shows the line of second (first) order 
deconfinement transitions. The facts that unlike the A = case, the 1st order line does not 
have an endpoint here, and the line coming from the large (3f side is clearly not a decon- 
fining line before the line coming from the Wilson axis meets it at (3f ~ 2, give credence to 
the hypothesis that the line in the large f3 v region is a coincident bulk and deconfinement 
transition line and, as we will see later, allows one to actually see two separate transitions 
on the same lattice. In what follows, we discuss first the 1st order transition line and then 
the deconfinement transition line in some more detail. 



A. The 1st Order Transition Line 



The first order transition points, shown in Fig. are listed in Table | along with the 
discontinuities in various observables and the values of at the transition point on the 

low (3f (P v ) side at fixed (3 V {(3f). In the f3 v > (3f region, the transition is associated with large 
discontinuities in P and (\Lf\). Also (\Lf\) ~ till the transition point, indicating that the 
transition line signals a deconfinement transition. The transition line in the large (3 V region 
is very similar to the A = case. This is expected, since in the large j3 v limit, the monopoles 
are automatically suppressed(see Fig. [1]) and the additional monopole suppression term 
does not have much effect. As j3 v — > 00, the dominant contribution to the path integral 
comes from those configurations which have a p .TrfU p = 2. Now the plaquette variables a p 
can be integrated out, constraining the gauge variables Ui to take values only in the center 
group Z 2 : o p = |Tr/[/ p = I\i e g p cri where 07 are Z 2 variables defined on the links and the 
Ui are frozen to the values o\. Then a c = 1 for all cubes c in this limit, and therefore the 
additional monopole suppression term does not have any effect. The action reduces to that 
of the Z 2 gauge theory, which has a well known first order transition at (3f ~ 0.44 with a 
discontinuity in the average plaquette. As j3 v is reduced from 00, the transition line starts 
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shifting from that of the A = case. The transition line shown in Fig. [| can be reproduced 
reasonably well up to /3 V ~ 2 by taking into account only configurations where U\ can have 
small fluctuations around 07: 

Ui = o\U[ where JJ[ — 1 + %gA\ . 
Perturbatively integrating over the fluctuations yield, in the leading order, 

Po= i< CTplw) = 1 --L + 0( _l). 

Up to this order, the integral over Ai only produces a renormalization of f3f, and the tran- 
sition point changes to 

^J^y (6) 

V 4/3 J 

The prediction of Eq. (||), shown in Fig. 0, matches the Monte Carlo results quite well for 

8 > 2 

In the large f3f region, the transition is associated with a large discontinuity in P a , as can 
be seen from Table |. Also from the (\Lf\) value at the lower side of the transition point one 
can see that both the sides of the transition are in deconfined state. The transition line in 
this region can be understood from the known bulk transition in the Z2 gauge-Higgs theory. 
In the (3f — > 00 limit, the gauge variables are frozen : Tt/U p = 2 Vp. The action in Eq. ([5]) 
then reduces to 

5 = A,E( 1 - ff p) + A E( 1 -^) , (7) 

p c 

after dropping an irrelevant constant. It describes a system with only Z2 degrees of freedom 
<jp. Under a duality transformation JI3| , one can rewrite it as the action for a Z 2 gauge-Higgs 
system (modulo irrelevant constants): 

S = -Pv E II 7i - ^ E s ili,i+» s i+» ■ ( 8 ) 
p ledp i,n 

Here s and 7 are Z2 variables residing on the sites and links of the dual lattice respectively, 
and 

^} = ilncoth{A. (9) 



The Z2 gauge-Higgs system is known to have a nontrivial phase diagram ||14|| ) leading to a 
nontrivial structure in the (3 V - A plane at /3f = 00 for the mixed Villain action. For A = 1, 
a first order transition at j3 v ~ 0.44 is predicted for 0f = 00. For large but finite (3f, the 
gauge variables fluctuate around the frozen value and yield the leading order estimate for 
the transition point, 

A - tt^Vy ■ ( 10 ) 

4/3/ ) 
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This leading order estimate is seen in Fig. || to work quite well up to 0f ~ 2. It is also 
consistent with Table |, which shows the transition to have a nearly continuous P and an 
almost constant discontinuity in P a in the region of very large f3f. 

While the above effective descriptions of the bulk transitions for f3f — > oo and f3 v — > oo 
were different, they are probably related by a duality transformation. In the A — > oo limit, 
the theory is self-dual ||15|| . In this limit a c is constrained to be 1, which can be solved as 
°v = T\i&p a u where o~i are Z 2 variables defined on links. Now a transformation U\ — > U{Oi 
interchanges the role of the fundamental and adjoint terms, leaving the action invariant. 
Under this duality transformation, flf <-> (3 V and P <-> P a . Even for A = 1, the constraint 
cr c = 1 is approximately true, and the self-duality is approximately obeyed. The locations 
of first order transition points in Fig. |] and the discontinuities in the plaquette variables 
P and P a presented in Table ^ display such a symmetry around the (3f = (3 V line. Taken 
together with the good agreement of Eqs. (|6]) and (|ToD with the Monte Carlo results in Fig. 
0, it leads one to the identification of the entire first order transition line traced out by the 
data points as a bulk transition line. Since the order parameter for the deconfinement phase 
transition, (\Lf\), becomes nonzero on the same line for f3 v ^ 0.74, one has a coincident 
deconfinement transition line as well. To investigate the origin of the deconfining nature of 
the first order line for f3 v > 0.74, we next turn to an exploration of the deconfinement line 
starting from the small (3f side. 



B. The 2nd Order Transition Lines 

For (3 V = 0, the chemical potential term decouples from the SU(2) part and we have the 
usual Ising-like second order deconfinement transition on this axis, at j3f ~ 2.3 for N r = 4 
lattices [|J. One expects, from continuity, a second order transition also for small V , as in 
all the examples mentioned in Sec. |. Hysteresis runs at (3 V = 0.3, 0.5 and 0.7 did indeed 
indicate a continuous deconfinement transition. For a more quantitative study, we studied 
the |L/ (-susceptibility, X\L f \, on N% x 4 lattices at each j3 v with N a = 8, 12 and 16. For 
each lattice, a low statistics run was made at the transition point estimated from hysteresis 
runs, and the peak of the susceptibility estimated by extrapolating to nearby couplings using 
Ferrenberg - Swendsen methods ||16|| . If the peak was not too far from the input /?/, a longer 
run was made to generate 10 5 configurations. Otherwise the procedure was repeated at the 
fresh peak location. The susceptibility curves for the three (3 V values are shown in Fig. || 
The critical exponent uj at each f3 v , obtained from a linear fit to In {x\L f \)max = &> In N a , is 
shown in Table |I| and is seen to be in good agreement with the (3 V = exponent. The average 
plaquette P from these runs was smooth everywhere, and the corresponding susceptibility 
peaks did not sharpen with N a at all, indicating a lack of bulk transition at these points. 
The transition points, shown in Fig. ||] by triangles, are therefore pure finite temperature 
transitions. 

On increasing j3 v further, the transition was found to change its behavior. Already at f3 v 
= 0.74, the rise of Polyakov loop from zero was associated with a discontinuity, indicating a 
first order deconfinement transition. The first order nature of the transition at this coupling 
was ascertained by a finite size scaling study. The distribution of the plaquette and Polyakov 
loop variables at the transition point are shown in Fig. |j. They show the presence of 
metastable states, and with increase of spatial lattice size, the two - peak structure was seen 
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to sharpen further without any visible movement in the peak positions, indicating a first 
order transition in the infinite volume limit. Fig. ^ also clearly shows that the transition is 
a deconfining one, since the Lf distribution at the lower state is peaked about zero. 

On increasing f3 v further, the same behavior continued : the deconfinement transition 
line, characterized by rise of from ~ to a nonzero value, was associated with a 

discontinuous jump. Also the average plaquette displayed a discontinuity at the same point. 
In Fig. ^ the complete deconfinement line is shown with the dotted part indicating second 
order and the solid part indicating first order transitions. This line is very much reminiscent 
of the deconfinement transition line of Refs. and is suggestive of a merger of the 

deconfinement transition line with the bulk transition line. 

A possible litmus test of the coincidence scenario is to see two separate transitions on 
the same lattice, especially in the vicinity where the two transition lines meet and remain 
coincident thereafter. From the transition points in Table | and Fig. ^| one can see that at 
(3 V ~ 0.7, one may see two separate transitions for N T = 4 lattices. Such an expectation was 
borne out by an explicit study at (3 V = 0.7 on an 8 3 x 4 lattice, the results of which are shown 
in Figs. |5]a and |^b. The former, a hysteresis run performed from the low /?/ end, shows 
that all observables, P, P a and (\Lf\) show a jump at around j3f ~ 2.21 whereas the latter 
shows the order parameter (\Lf\) along with its susceptibility, obtained from the longer run 
mentioned above. It clearly shows a second order deconfinement transition taking place first 
at Pf ~ 2.1, followed by a bulk phase transition later at (3f ~ 2.2. 

We have carried out a similar exercise for N T = 6 lattices as well, with essentially similar 
results. The deconfinement trajectory now starts from the Wilson axis at f3f m 2.43 |)17| . At 
(3 V = 0.3 and 0.5 we get a continuous transition, with the transition point shifted to a slightly 
higher (3f (see Fig. |), as expected for a physical transition. The susceptibility curves for 
N a = 12, 14 and 18 lattices, shown in Fig. |B|, and the critical exponent obtained from 
them, shown in Table |J, indicate an Ising - like second order transition at these points. On 
increasing j3 v further, the transition line hits the 1st order line and turns first order, but now 
at a slightly smaller value of (3 V . Already at f3 v = 0.7 one sees a first order transition, as can 
be seen form the plaquette and Polyakov loop distributions shown in Fig. 0. One thus needs 
to choose a slightly smaller value of (3 V for N T = Q in order to see again two transitions on 
the same lattice but one does see them. The A = 1 simulations for the mixed Villain action 
thus lend a strong credibility to the hypothesis that the deconfinement transition line for 
possibly a large range of N T merges with the bulk transition line. Since the latter branches 
out in this case and exhibits no end point, the merger is easy to observe numerically: the 
small f} v region has only a deconfinement transition line and the large (3f region has only a 
bulk transition line, while they seem to be coincident in the (3 V ^ 0.7 region. 

While the Z2 - monopole suppression thus provides convincing indications that the un- 
physical bulk transition line is coincident with the physical deconfinement transition line for 
a range of N T , and thus provides a reasonable explanation for the latter turning first order 
in the f3f - j3 v plane, it is insufficient to remove all the bulk transitions from the region of 
interest and to show a second order transition similar to the Wilson action case for large (3 V 
values. We therefore turn our attention in the next section to the other Z2 objects causing 
the remaining bulk transition : the electric vortices [0. 
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IV. SUPPRESSION OF ELECTRIC VORTICES 
AND MAGNETIC MONOPOLES 

In order to suppress the electric vortices in addition to the magnetic monopoles, we 
consider the action 

S(U,a) = f3 f J2(l- 1 -Tr f U p S ) 

+&e( i -^-' i w) 

+A£(i-^)+7E(i-^), (ii) 

c I 

with a e defined in Eq. (£|). The last term suppresses the Z 2 electric current loops. It should 
be irrelevant in the continuum limit, leaving the relation between the continuum coupling 
and f3f and (3 V unchanged. For sufficiently large 7, one expects that the bulk transition line 
caused by the condensation of electric loops will be suppressed. One can reduce the number 
of input parameters by sending A — > 00 and using the the explicit self-dual form of the above 
action. The monopole term is then absent and the plaquette variables a p are replaced by 
products over corresponding Z 2 - link variables. This form has the advantage of being more 
tractable numerically. We employed this self-dual form for our investigations at large 7. 

We followed the same procedure as in the previous section to monitor (\Lf\) and the finite 
size scaling behavior of its susceptibility for studying the finite temperature transition. P 
and P a were utilized to look for the bulk transitions. We used a heat-bath algorithm for 
both the gauge and Z 2 variables. One sweep consisted of updating all the gauge variables 
of the lattice, followed by updating of all the Z 2 variables. We made a number of hysteresis 
runs of 15000 iterations per point at different (3 V values on a 8 3 x 4 lattice. No discontinuous 
transition was found anywhere in the f3 v — (3f plane. The plaquette variables P and P a were 
smooth everywhere, indicating that the additional term has succeeded in eliminating all bulk 
transitions. The deconfinement transition points were estimated from the rise in (\Lf\). In 
order to look for the expected shift in the deconfinement transition point, we repeated the 
exercise on a 12 3 x 6 lattice. Our investigations revealed lines of deconfinement transition 
points, consistent with (see Fig. 

p f + A « , (12) 

where 0^ is the deconfinement transition point for the Wilson action (« 2.30 for iV T =4 |§ 
and ~ 2.43 for N T = 6 []T7| lattices). Rapid but continuous rises of (\Lf\), associated with 
diverging susceptibilities with increasing spatial volume, were found at the deconfinement 
transitions at all the f3 v values investigated. These properties are similar to the known 
features of the deconfinement transition for the Wilson action for SU(2). 

The critical exponent uo was obtained from a finite size scaling analysis of the In- 
susceptibility on iV 3 x 4 lattices with N a = 8, 12 and 14 at f3 v — 1 and 2. Long simulation 
runs of 10 5 thermalized configurations and the usual spectral density methods were used to 
obtain the peak height at the transition point for each lattice. The susceptibility peaks are 
shown in Fig. |9|. For iV T = 6, we used N a = 12, 14 and 18 and investigated the f3 v values 0.5 
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and 1.5. The results are shown in Fig. |10[ Values of the critical exponent uj obtained from 
these peaks are listed in Table III, and are in good agreement with the value for the Wilson 
action for SU(2). This behavior is, of course, what one expects from universality arguments. 
Moreover, the phase transition line, Eq. (|12D , is also consistent with the expected continuum 
limit behavior of this action. 

Suppressing the artifacts which cause the (first order) bulk transitions in the mixed 
Villain action, namely, the Zi magnetic monopoles and vortices defined above, thus yield 
results consistent with expectations from universality, both for the order of the transition as 
well as the shift of the transition point with iV r . Perturbatively, one expects the contribution 
of these topological objects to be zero in the continuum limit. It is, however, unclear 
whether this is indeed so nonperturbatively as well. In fact, the earlier results for the 
unsuppressed action (A = 7 = 0) suggest that very large N T ( 3> 16 ) will be needed to see 
the deconfinement transition separate out from the bulk transition in the f3f - (3 V plane. 



V. SUMMARY AND CONCLUSIONS 

Earlier numerical simulations on asymmetric N% x N T lattices using mixed actions, de- 
fined by Eq. P) or (^), yielded apparently paradoxical results. The deconfinement phase 
transition in these cases was found to become first order for a large range of the adjoint 
coupling, although one naively expects these actions to be in the same universality class as 
the Wilson action, which is known to have a second order deconfinement transition. On the 
other hand, the transition exhibited very little shift with the temporal lattice size N T and 
the average action was also discontinuous at the transition. These results could be consistent 
with universality if a curious interplay between the physical deconfinement transition and 
the bulk transitions, which are merely lattice artifacts, existed for some unknown reasons. 
The absence of any significant variations in the locations of the transition or any change in 
its nature for large f3 v made other options of saving universality rather implausible. 

Since variations with N T , a defining characteristic of the deconfinement phase transition 
which is not supposed to affect the bulk transition, failed to resolve the paradoxes, it seemed 
natural to seek a mechanism to shift or eliminate the bulk transition(s), leaving the univer- 
sality class intact. Taking a clue from our earlier studies for the SO(3) gauge theory, where 
a similar mixture of bulk and deconfinement effects is seen, we attempted in this paper 
elimination of the bulk transition(s) by suppressing certain topological objects as a means 
to study the deconfinement transition in the mixed theory. 

First, we considered the mixed Villain action with an addition of a chemical potential 
for the Z2 - magnetic monopoles and studied the phase diagram on asymmetric lattices 
in the f3f - f3 v plane of couplings. The monopoles get suppressed with increasing chemical 
potential A. Numerical simulations for A = 1 showed an interesting phase diagram which was 
quite different from that of the original theory. Nevertheless, it still shared the paradoxes 
mentioned above. We could, however, demonstrate that i) a second order deconfinement 
transition line emanating from the (3f axis meets the bulk transition line at a finite (3 V , ii) 
the bulk transition line has no end point, and iii) the change of the order of the deconfinement 
phase transition occurs as the two lines merge. Moreover, we were able to show the presence 
of two phase transitions on the same finite lattice in the vicinity of the point of merger : 
a second order deconfinement phase transition, indicated by a continuous sharp rise in the 
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order parameter, followed by a discontinuous bulk transition, where all observables exhibited 
an abrupt change, but both sides of the transition were in a deconfmed state. 

A further suppression of the Z 2 electric vortices enabled us to get rid of the bulk transi- 
tions completely, as there was no trace of any discontinuity anywhere in the coupling space 
and the average plaquette variables were smooth everywhere. A study of deconfinement 
transition in this theory revealed a line of second order deconfinement transitions, obeying 
/3f + /3 v ~ (3^1 where 0^ ls the corresponding transition point for the Wilson action. There- 
fore the removal of the bulk transition lines by suppressing the topological objects causing 
those transitions led to results fully consistent with universality. Since the terms added to 
the action in the process do not contribute in the naive continuum limit, one can formally 
attribute the anomalous behavior of the deconfinement transition lines for the action in Eq. 
(0) to the presence of bulk transitions. Of course, it is natural to expect that the same thing 
happens for the action in Eq. ([!]), too. While the vexed possibility of violation of universal- 
ity in these cases is thus finally eliminated, some nagging questions still remain. Foremost 
amongst them is how rapidly the topological objects discussed here vanish from the lattice 
for the original theories as one approaches the continuum limit. Simply put, even for the 
original theories one now expects for sure the separation of the two lines of transitions for 
some N T . However, it is also clear that it must be at an almost astronomically large value. 
One knows that this is indeed the case for spin models [fTB|| , where one needs to work at 
a billion times larger correlation lengths to see universal behavior in the adjoint direction. 
One does not expect it to be so, however, for a continuous gauge symmetric case as ours. 
Another interesting and related question is about the behavior of the order parameter at the 
bulk transition when the two transitions do separate : it will have to be continuous and zero 
through the bulk transition and it will also be exponentially small after the deconfinement 
phase transition for very large N T due to the divergences associated with the point source 
it represents. 
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FIG. 1. The phase diagram for the action (y), showing the first order bulk phase transition 
lines (from Ref. ||). The bulk phases are characterized by values of E and M, as shown. 
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FIG. 2. The phase diagram for the action (|5|) with A = 1 on an 8 3 x 4 lattice. The filled circles 
show first order transition points. The dashed lines show the estimates from Eqs. (0) and (|To|). 
The triangles in the low f3 v region show the locations of Ising-like second order deconfinement phase 
transitions on N T = 4 and 6 lattices respectively. The dotted and solid lines show the second and 
first order deconfinement transition lines. 
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FIG. 3. X|L/| as a function of 0f for iV T = 4 lattices with iV CT = 8, 12 and 16 for the action 
(gj) at /3„ = 0.7, 0.5 and 0.3. The values obtained directly from the long simulation runs are also 
shown with error bars. The lines were obtained using Ferrenberg-Swendsen extrapolation. 
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FIG. 4. Probability distributions of Lf and the plaquette variable at (3 V = 0.74, /3f = 2, for N T 

= 4 lattices with N a = 8, 12 and 16. Two peak structures are visible that sharpen with increase 
in lattice size, indicating a first order transition. The Lf peaks clearly indicate that the transition 
is a deconfining one. 
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FIG. 5. (a) Results of a hysteresis run at j3 v = 0.7. The diamonds, triangles and circles denote 
the variables P, P a and (\Lf\) respectively. (\Lf\) is seen to rise at Pf ~ 2.1, and a discontinuous 
transition is seen at Pf ~ 2.2. (b) and its susceptibility as a function of Pf at f3 v = 0.7, 

showing clearly a second order deconfinement transition at /?/ ~ 2.1 and a first order transition at 
f3 f ~ 2.2. 
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FIG. 6. X\Lf\ as a function of /3j for iV T = 6 lattices with N a = 12, 14 and 18, for the action 
(gj) at (3 V = 0.5 and 0.3. The values obtained directly from the long simulation runs are also shown 
with error bars. The lines were obtained using Ferrenberg-Swendsen extrapolation. 
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FIG. 7. Probability distributions of \Lf \ and the plaquette variable at j3 v = 0.7, j3f = 2.12 
for N T = 6 lattices with N a = 12, 14 and 18, showing metastable states indicating a first order 
transition. 
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FIG. 8. Locations of the deconfinement transition points in the (/3f , (3 V ) plane. The circles and 
diamonds show the simulation points for N T = 4 and 6 lattices, respectively. The lines correspond 
to Eq. (0). 
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FIG. 9. Polyakov loop susceptibilities for N% x 4 lattices for N a = 8, 12 and 16, for (3 V = 1 and 
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FIG. 10. Polyakov loop susceptibilities for N^. x 6 lattices for N a = 12, 14, 18, for (3 V = 0.5 and 
1.5. The errors are shown for the simulation points. 
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TABLES 

TABLE I. Transition points and discontinuities for the first order transitions in the monopolc 
suppressed mixed action. 



Pvc 


Pfc 


AP 


AP a 


A(|L/|) 


\Lf\- 


5.0 


0.52(2) 


0.36(2) 


0.007(3) 


0.60(4) 


0.04(3) 


3.0 


0.60(2) 


0.33(2) 


0.022(4) 


0.46(4) 


0.03(2) 


2.0 


0.76(2) 


0.26(1) 


0.08(1) 


0.38(6) 


0.07(5) 


1.0 


1.60(4) 


0.17(4) 


0.22(5) 


0.35(6) 


0.06(4) 


0.74(2) 


2.0 


0.11(1) 


0.32(1) 


0.33(6) 


0.05(4) 


0.59(2) 


3.0 


0.023(4) 


0.33(2) 


0.04(5) 


0.51(4) 


0.52(2) 


5.0 


0.007(2) 


0.36(2) 


0.01(6) 


0.69(4) 



TABLE II. Values of the critical index to for the monopole suppressed action. 





N T = A 


N T = 6 




0.3 


0.5 


0.7 


0.3 


0.5 


UJ 


1.92 ±0.05 


1.87 ±0.05 


1.91 ±0.02 


1.82 ±0.08 


1.89 ±0.07 



TABLE III. Values of the critical index to for the action with both monopoles and vortices 
suppressed. 





N T = 4 


N T = 6 


Pv 


1.0 


2.0 


0.5 


1.5 


UJ 


1.88 ±0.05 


1.93 ±0.05 


1.84 ±0.06 


1.86 ±0.07 
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